Hirota equation as an example of integrable symplectic map[| 
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Abstract. The hamiltonian formalism is developed for the Sine-Gordon model on the space- 
time light-like lattice, first introduced by Hirota. The evolution operator is explicitely con- 
structed in the quantum variant of the model, the integrability of corresponding classical 
finite-dimensional system is established. 

Introduction 

We consider a system of difference equations on a plane lattice shown below. 




The variables, attached to the vertices of lattice, are denoted by w. The equations connect the values 
of iy's belonging to a common elementary square 



UK 



-s^f (1) 



as follows 

/K 

The equation allows to get w N if three other tu's are given and so defines an evolution in the vertical 
direction (from south to north). 

If the function f(w) is invertible, so that the equation 

w = f(w) (2) 
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has a unique solution 

w = g(w) (3) 
then the equation (||) can be rewritten in the form 



9(w s ) 
: 5(* N ) 



(4) 



and defines the evolution from east to west. In what follows we shall consider only move south-north, 
defining the time axis as vertical. Then the horisontal saw 



can be considered as a "line" of initial data. 

It is evident that the initial data, given on a portion of that saw, define a solution on the part of a 
lattice bounded by the punctured lines: 



In this sense we can speak of hyperbolicity of our system, the lines being lightlike and defining the light 
cone. 

Alternatively we can use a kind of dual variables ip to define the same dynamical system. These 
variables are attached to the vertices of a "dual" lattice 





The connection of ?/>'s and w^s is given by 

V. 

w = ^ (5) 

and 

iP = l[w (6) 

where the product is taken over all tu's on the same horisontal line to the east of ip. The second formula 





needs clarification of the boundary conditions; this will be done in the main text for the lattice periodic 
in the horisontal (spacial) direction. 

In terms of ^'s the equations take the form 



which is also eqivalent to 



^ N =V s /(^) (7) 



(see the last picture). 

In Section 1 we present a symplectic structure on the space of initial data which is conserved by 
evolution and construct the generator of this evolution. We shall do the latter in a quantum version 
of our system where this generator is given by an unitary operator. Then in Section 2 for a concrete 
choice of the function / 

tt \ 1 + k2w rn \ 

/H = -jt — (9 

where k 2 is a real positive parameter we shall show that our symplectic map is completely integrable. 
Namely, we shall present the set of independent commuting integrals of motion, as many as the number 
of degrees of freedom. Finally in Section 3 we shall show the equivalence of our system to that introduced 
by Hirota |l| as the one approximating the famous sine-Gordon equation. Some general comments are 
collected in Conclusion. Some results of the present paper appeared already in ||. 

1 Symplectic map 

We shall consider our system on a spacially periodic lattice of even length 2M. The variables wona 
saw of initial data will be numbered as shown on the picture: 



W2n+1 W2n 




W2n-2 



We introduce the following Poisson structure 

{w n ~i,w n } = 2jw n -iw n n = 2,3, ... ,2M 

{W2M,W\} = 2JW2MWI (10) 

{w m , w n } = 1 < \m — 7i | < 2M — 1. 

Here 7 is a positive parameter playing the role of a coupling constant. The Jacoby identity is evident 
as the Poisson matrix becomes constant after passing to logarithms of tu's. 

In applications one can encounter different variants for the set of values of the variables w. The 
most evident are the following ones: 

1) w real and positive; 

2) w complex and unimodular, \w\ = 1. 

Poisson brackets are compatible with both cases; there are also other possibilities but in this paper we 
confine ourselves to the second (compact) form of w. It is thus necessary to require that the function / 
maps the unit circle on itself. 

The Poisson structure is degenerate; the functions 

C\ = l\ w odd ^-q 



are commuting with all 10's. The submanifolds with fixed values of C\ and C2 realize the symplectic 
leaves of our manifold. The number of degrees of freedom is thus M — 1. 
The map 

U : w — ► w 

where w's are the values of the solution of the equations ([l]) on a saw S shifted by time unit from the 
initial saw S 



W 2n +1 W 2n -l 



s 



W2r 



W 2n -2 



W 2n +1 W 2n -1 

\/\/\/\ 

W 2 : 



s 



W2n-2 



preserves the Poisson brackets and the values of invariants C\ and C2. Thus IA defines a symplectic map 
(canonical transformation). It can be described by means of corresponding generating function. We 
shall not follow this route here and instead shall describe the unitary evolution operator U corresponding 
to this map in the quantized variant of our dynamical problem. The reason for this is, that quantum 
mechanics being more algebraic, the corresponding formulae are just more simple. One can then get 
the classical variant by the appropriate limit. 

After quantization the variables w become unitary operators with the commutation relations 

w n -iw n = q 2 w n w n - 1 n = 2, 3, . . . , 2M 

W2MW\ = q 2 w 1 w 2 M (12) 
w m w n = w n w m 1 < \m — n\ < 2M — 1. 



where 



h being the Planck constant. This commutation relation passes into the Poisson relation (|lG|) with 
h — > if we identify as usually |[ , ] = { , }. 
The evolution operator is given by 

U = JJ r(w even ) Y[ r(w odd ) (13) 

where the function r(w) is defined by the function f(w) through the functional equation 

= /(-)■ I") 

r(q l w) 

There is no ordering of factors problem here as only one variable enters the functional equation while 
in (13) we took into account that all even w's commute among themselves, the same for odd ones. 

Application of U to w's is rather trivial due to the locality of their commutation relations. For even 
w's we get 

W2n = U~ 1 w 2n ~U = r~ x (w 2n -i (w 2n +i ) w 2n r (w 2n +i )r (10271-1 ) (15) 
which after using the functional equation (|l4|) gives 

W2n = f(qW2n+l)'W2n(f(qW2n-l)y 1 - (16) 

For odd w's the calculation consists of several steps but it is as trivial as above. We get 

Vhn-l = U~ 1 W 2 n-lU = f {qW2n)W2n-l(f {qW2n-2)Y l ■ (17) 



So the iteration of U leads to the quantum form of the equations of motion 

w N = fiqw^Wsifiqw^y 1 (18) 

which coincides with (|l]) in the classical limit q — > 1. 

Consider now the alternative picture in the variables ip. The set of variables ijj n 

1p2n 1p2n-2 




will be now quasiperiodic 



V"2n+1 1p2n-l 



(19) 

^271-1+2^" = ^2^2n-l+2M 



and 2M + 2 independent variables C\, C2, 1(^1, tp2, ■ ■ ■ , '4 ) 2M (M + 1 degrees of freedom) have the following 
system of commutation relations: 

V'mV'n = q^n^m, if n — m odd and < n — m < 2M 
ipoddCi = q 2 Clip odd 

i>evenC2 = q 2 C 2 1peven (20) 
VPoddAodd] = [ip 'even A 'even] = 0. 

Introduce the evolution by means of the same operator U with 

W n = ~ • (21) 



Of course these ui's define system ( J12D with C\ and C2 as monodromies (|TT|) . 

The derivation of equations of motion for ^'s is even simpler than that for w's because each ip does 
not commute with only one w 

ip n w n = q 2 w n i\) n . (22) 



We get 

^ = ^ s f{q- 1 ^) rill) 



which passes to (0) in the classical limit q — > 1. 



2 Integrability 

We consider now the particular case with the function f(w) given by 

1 + K 2 W 

fW = 2 1 • 24 
K + W 

Here k 2 is a fixed positive parameter. We thus have |/[ = 1 if |u>| = 1. The inverse function g is given 
by the same formula after reflection k 2 — > —k 2 . 

Of course this function has a very special character and it is not clear a priori why it is interesting. 
Its exclusiveness is due to the following property ||] : the solution of the functional equation 

r(X,gw) 1 + Xw 



r(X,q 1 w) \ + w 



which we use in the definition of the evolution operator ([13]) 

u = II r(K 2 ,w even ) JJ r(K 2 , Wodd) 
satisfies the Yang-Baxter equation in the following form: let u and v be a Weyl pair 



uv = q vu, 



(26) 



(27) 



then 



r(A, u)r(\/i, v)r({J,, u) = r(fi, v)r(\fj,, u)r(X, v). (28) 

This equation allows to interprete r as a fundamental L-operator of quantum chain 0j and use the 
machinery of QISM (see e.g. ||). In particular we can introduce a family of matrix L-operators with 
operator matrix elements 

/ 



1 i 

2 .1,2 



Due to (^) it is connected with r by the commutation relation H 

r(A «, n )i n+ i(A)L ft _iOx) = V§M L n-§( A )K^>n) 

where presumably w n = f" +1 . 

One can recognize in L a Lax operator for the Hirota equation. Indeed, the equation 
interpreted as a zero curvature condition for the elementary square of our lattice: 

/ \ 

L w< _ n (k x A) L n ^ e (/«A) 

\ / 

L w< _ s (kA) L s ^ e (k *A) 



(29) 



(30) 



could be 



where 



L w< _ n (k 1 A)L n ^ e (kA) = L w+ _ s (kA)L s ^e(k X A) 
ip A 2 tpi XipA 2 ipn 2 ] 



i i 

, 2 „/. 2 



1 _1 

2 »/, 2 



(31) 



(32) 



V AV« V'i^B 2 / 



It's basically the same zero curvature representation as the one invented by Hirota [GJ but this time it 
is truly quantum. 

Multiplying the -L-operators along the initial data saw and adding the additional diagonal factor to 
take care of the quasiperiodicity of ?/>'s (see [Q]) we get the monodromy matrix 



T(A) = L 2M _i(/eA)L 2M _3 (k _1 A) . . . L|(kA)Li (k _1 A)C£ 3 
with commuting traces t(A) = trT(A) 

[i(A),t(/i)]=0. 



(33) 



(34) 



The last statement is not quite obvious because the L-operator does not seem to fit into the standard 
ultralocal scheme of QISM. It's somewhat similar to the L-operator of the quantum Volterra model || 



but involves highly nonlocal variables ip. The ultralocality is however easy to recover by glueing our 
L-operators in pairs 



L„(A) 



L 2n _i(K\)L 
( 



X 



t(X) 



A VL-1 + A ™2n-1 

1 I _I _I 

V ^2n^l( K_1 Wn-l + Kw 2n-l)^2n-l 

tr(L M (A)L M -i(A)...L 1 (A)C 2 ff3 ). 



l 

2 

2n-l 



(35) 



Looking at the commutation relations for the operators involved 

1p2n-lW2n-l = q 2 W 2n -1^2n-^ 



bp2m-l, V)2n-l] 



if m 7^ n 

[w2m-l,W2n-l] 



(36) 



0. 



we notice that L n is not only ultralocal (i.e. matrix elements of the operators L n with different numbers 
n commute) but it also resembles the L-operator of the general lattice XX Z model (compare also with 
the lattice sine-Gordon L-operator [@]) and indeed one can easily verify that 



fl(£)L n (A) ® L n (fi) = L n (/x) ® L n (A)i?(^) 



(37) 



where R(X) is the usual 4x4 trigonometric i?-matrix Q. This makes the commutativity of t(X) almost 
evident while the commutation relation (p0|) (or the zero curvature equation (|3l|)) allows to establish 
easily that t(A) is conserved by our dynamics 



U- x t(A)U = t(X). 
Direct inspection shows that t(X) is a polynomial of A 2 of degree M 



M 



t(X) = 4A 

fc=0 



2 A- 



(38) 



(39) 



and thus M + 1 coefficients Ik constitute the number of commuting integrals of motion required for the 
complete integrability. It remains to establish their functional independence; this can be done by usual 
methods. 

In classical limit these integrals remain nontrivial and their existence proves the integrability of our 
system. 



3 Hirota equation 

The equations considered above do not coincide exactly with the Hirota equations. The correspondence 
can be achieved by means of the following change of variables: on each second diagonal line SW-NE of 
the (dual) lattice 



we are to make change tp — > tp .If the function f(w) is "odd" 

f(w- 1 ) = r\w) (40) 

(the function f(w) from Section 2 is "odd" while the function r(A, w)) is "even": r(X 1 w^ 1 ) = r(X,w)) 
then the equations of motion in the new variables (call them x) take the form 

x N = x^/OT^wXe)- ( 41 ) 

Instead of w 7 s we now introduce it's 

u = X w Xe- ( 42 ) 

They are connected exactly in the same way as ^'s and namely, u coincides with either w or w^ 1 
on alternating diagonals of their lattice. The equations of motion in terms of u's read 

u N = f(qu w )u- l f(q- l u E ). (43) 

Of course on a periodic lattice all these changes are possible only if the period is even. 

The commutation relations for the initial data and the evolution operator now look as follows: 

U2n±lU2n = q 2 U 2 nU2n±l U = 1, 2, . . . , M - 1 

U\U2M = q 2 U 2 MUl ^ 
U2M-1U2M = Q 2 U2MU2M-1 

u m u n = u n u m 1 < \m — n\ < 2M — 1. 

U = U Y[r(K 2 ,u even ) Y[r(K 2 ,u odd ) (45) 
where Uo is a reflection of all operators u n 

XJu 1 u n XJ = u- 1 . (46) 



The commutation relations ( |44| ) are compatible with this inversion which ensures the existence of Uo- 
The conservation laws I n are even function of u>'s and so are commuting with Uo- 

It is instructive to remind how equation (|4l]) turns into the sine-Gordon equation in the continuous 
limit. Introducing variables if 

X = e l * (47) 
we come to the Hirota form of the equation ( |4l| ) 

sin K<p w + tp E - <p N - ip s ) = kT 2 sin ±(<p w + tp E + ip N + <p s ). (48) 

Suppose that when the lattice spacing A goes to zero the variables 92 turn into the smooth function 
ip(t, x) and thus 

A 2 

Vw + Vv-Vk-Vs -^^(dlv- 9 ??)- ( 49 ) 

If now 

1 ™2 

(50) 



A 2 k 2 4 

then for (p(t, x) we obtain exactly the sine-Gordon equation 

2 

a t V-4V+^-sm2^ = 0. (51) 

Of course, this limiting procedure is valid only in classical case. The quantum case requires a more 
delicate balance between the dimensional parameter A and dimensionless k 2 , which accounts for the 
mass renormalization. 



Conclusion 



Depending on the interests of a reader, one can look upon this paper from different points of view. For 
a specialist in dynamical systems the Hirota model gives a family of examples of symplectic maps of 
a finite-dimensional phase space. In particular, as was explicitely mentioned in || for one degree of 
freedom we get an integrable deformation 



x n+ i - 2x n + x n -i = -2i In 



k 2 + e" 



k, 2 + e tXn 
of the so-called "standard map" 

x n+1 - 2x n + x n _i = -4k sin x n . 

For related results on symplectic maps see also || 10 1. 

For a specialist in quantum field theory the quantized Hirota model defines a lattice regularization 
of a relativistic local field-theoretical model, namely sine-Gordon model. This gives another example of 
the relation between alternating chain and relativistic massive models observed earlier in pi] , |i~2fj . We 
give a reasonably explicit construction of the evolution operator U = e~ lH ^ 1 where A t is lattice spacing 
in time direction and H is the hamiltonian. Inclusion of U into algebra generated by local Lax operator 
allows to introduce the algebraic Bethe-Ansatz technique for the diagonalization of U. Here we enter 
the well developed domain (see Q and the references therein). 

A specialist in exactly solvable plane models of statistical mechanics should see here new way of 
writing and arranging of formulas known in the theory of Z n model. It is enough to stress that for q 
being a root of unity r(\,w) is exactly the Fateev-Zamolodchikov i?-matrix [[14|]. 

Finally, a reader interested in the infinite-dimensional algebraic structures is to see here fragments 
of a future theory of the algebra U q {sl{2)) with which we connect our hopes for the construction of 
correlation functions in integrable models of quantum field theory. This programm is being developed 
now by the Kyoto group [pj[|. 
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